Giant capacitance of a plane capacitor with a two-dimensional electron gas in a 

magnetic field 



Brian Skinner and B. I. Shklovskii 

Fine Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 55455, USA 

(Dated: December 10, 2012) 

If a clean two-dimensional electron gas (2DEG) with small concentration n comprises one electrode 
of a plane capacitor, the resulting capacitance C can be larger than the "geometric capacitance" C g 
determined by the physical separation d between electrodes. A recent paper^ argued that when the 
effective Bohr radius as of the 2DEG satisfies as <C d, one can achieve C ^> C g at low concentration 
nd 2 -C 1. Here we show that even for devices with as > d, including graphene, for which as is 
effectively infinite, one also arrives at C 2> C g at low electron concentration if there is a strong 
perpendicular magnetic field. 



I. INTRODUCTION 

In a standard parallel-plate capacitor, the capacitance 
per unit area C is equal to the "geometric capacitance" 
C g = e/And (in Gaussian units), where e is the dielectric 
constant of the medium separating the two plates and 
d is the distance between them. The relation C = C g 
is correct when both electrodes are made from a perfect 
metal, which by definition has an infinite electron den- 
sity of states. On the other hand, when one of the ca- 
pacitor electrodes has a finite density of states, there is 
an additional "quantum capacitance" contribution to the 
total capacitance that reflects the finite thermodynamic 
compressibility of the electronic charge. This quantum 
capacitance C q can be written C q — e 2 dn/d[j,, where e 
is the electron charge, n is the two-dimensional electron 
concentration on the electrode surface, and /z is the elec- 
tron chemical potential. Thus, the total capacitance of 
a capacitor containing one ideal metal electrode and one 
electrode with finite density of states can be written 



More generally, this capacitance is defined as C — 
edn/dV, where V is the electrochemical difference be- 
tween the two electrodes, maintained by some external 
voltage source. In the following discussion it is convenient 
to define the effective capacitor thickness d* = e/4irC, so 
that when the thermodynamic density of states (TDOS) 
dn/dfi is positive, the effective thickness d* of the capac- 
itor is larger than the physical thickness d. 

On the other hand, it has long been understood that 
for a low-density two-dimensional electron gas (2DEG), 
the TDOS can be negativ^^, owing to the emergence 
of strong positional correlations between electrons. This 
implies a negative quantum capacitance, or in other 
words, d* < d, as wa s first measured experimentally 
over two decades agcP®^l Recent experiments on ultra- 
thin capacitor devices have reported a very large negative 
quantum capacitance, so that d* jd is as small a d 20 * 21 * 0.6. 
These experiments, along with the recent development of 
nanometer-thick graphene capacitor device d 22 " 24 !, bring 
to the forefront an important and fundamental question: 
how small can d* be? 



In a recent series of papers^H, we showed that when 
the electron density is sufficiently low that nd 2 <C 1, the 
metal-2DEG capacitor can be thought of as a collection 
of electron-image charge dipoles (as shown schematically 
in Fig.Jl]), whose interaction acquires the form ~ e 2 d 2 /r 3 
due to the strong screening effect of the metal. As a re- 
sult of this screening, the capacitance becomes greatly 
enhanced above the geometric value at nd 2 <C 1 when- 
ever the effective Bohr radius as is much smallei^ than 
d. In the limit n <C aj^/d 4 , positional correlations be- 
tween electrons are lost due to the short-ranged nature 
of the dipole interaction, and the capacitance saturates at 
a finite value such that d* = /4 (for spin-unpolarized 
electrons) . 




FIG. 1. (Color online) Schematic picture of a capacitor made 
with a low-density 2DEG next to a metal electrode in the pres- 
ence of a transverse magnetic field. Electrons (filled circles) 
create image charges (open circles) in the metal electrode. 
The transverse magnetic field confines electrons to their low- 
est Landau level wavefunction with characteristic lateral size 
Ib, as shown schematically for the central electron by the red 
shaded area. 

One can expect, then, that for 2DEGs with large effec- 
tive Bohr radius > (I, there can be no enhancement 
of the capacitance above the geometric value. In par- 
ticular, capacitors made from graphene, for which as is 
effectively infinite due to the massless Dirac spectrum, 
seemingly should not produce d* < d. Such thinking is 
consistent with recent e xperi ments probing the quantum 
capacitance of graphen e 127 * 28 ^, which showed that the ca- 
pacitance is everywhere smaller than the geometric value 
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and, in the absence of disorder, tends to zero in the limit 

In this paper, however, we show that even for 2DEGs 
with large (or infinite) Bohr radius, an enormous en- 
hancement of the capacitance is possible when the ca- 
pacitor is placed in a strong perpendicular magnetic field. 
Such an applied field helps to confine electrons laterally 
and preserve their strong positional correlations even in 
the limitP^n — > 0, so that the capacitance diverges. This 
situation is illustrated schematically in Fig. [T] 

In this paper we focus our attention specifically on the 
lowest Landau level (LLL), where the filling factor v < 1, 
and we calculate d* as a function of v in the absence 
of disorder. Our primary result is illustrated in Fig. [2] 
which shows d* /d as a function of filling factor for differ- 
ent values of d/ls, where Is = y/hc/eB is the magnetic 
length. From Fig. [2] one can see that when the capacitor 
is thin enough that d/ls <C 1, the capacitance is greatly 
enhanced above the geometric value throughout the LLL. 
For large d/ls, this enhancement is large only at v <C 1 
and 1 — v <C 1. These results are valid for a capacitor 
made with either a conventional spin-polarized 2DEG or 
with a graphcnc layer where both spin and valley degen- 
eracies are lifted, due to an exact correspondence between 
the two systems in the LLL. We assume throughout this 
paper that we deal with such a nondegenerate LLL, so 
that the filling factor v = 2nnlg. 



LLL in cases where d/ls <C 1. Throughout this paper we 
work in the range of temperatures where the thermal en- 
ergy is small compared to the interaction energy per elec- 
tron, yet still large enough to smear out the singularities 
in energy related to the fractional quantum Hall effect. 
The effects of finite disorder are discussed in Sec. IIVI 
Finally, in Sec. [V] we discuss capacitors made from two 
parallel, low-density 2DEGs, where very similar physics 
also produces large capacitance enhancement. 



II. GENERAL NUMERICAL RESULT 

The general, thermodynamic definition of the differ- 
ential capacitance is C _1 = (l/e 2 )d 2 [nE(n)]/dn 2 , where 
E{n) is the total energy per electron. The effective thick- 
ness d* can thus be written 



rl 2 d 2 



Particle-hole symmetry in the LLlPl implies that E{v) 
must satisfy v[E{v) - E{\)] = (1 - u)[E(l - v) - E(l)}. 
As a consequence, d* in the LLL is symmetric about v — 
1 /2 (as can be seen in Fig. [2]) . Given the particle- hole 
symmetry requirement, FO suggested that the energy in 
the LLL can be written using the following power law 
expansion: 
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FIG. 2. The effective capacitor thickness d* , in units of d, 
plotted as a function of filling factor v in the LLL at different 
values of djls- 

In the remainder of this paper we first present our 
derivation of the capacitance (Sec. [Tl]), which is based 
on generalizing the Fano-Ortolani (FO) expression^ for 
the energy of a 2DEG in the LLL to account for screening 
by the gate electrode. We then discuss, in Sec. [TTT] the 
limit of large electron density, nd 2 1, at which the gate 
screening effect is unimportant, and we show that we can 
reproduce the traditional FO result. The opposite limit 
nd 2 <C 1 is also considered in detail, and we derive an 
analytical expression for d* that is valid throughout the 



vE(v) = E{iy + Y,<*hW(l-'')] k/2 

k=0 



(3) 



The condition E(0) = implies that ao = u\ = oti = 0. 

The energy E(l) can be calculated^ by integrating the 
electron-electron interaction law V(r) weighted with the 
pair distribution function 



g (r) = 1 - exp 



2Z 2 



(4) 



of the incompressible state at v = 1. For electrons with 
a nearby metallic gate, the interaction law is given by 



V{r) 



e 

sr 



:^r 2 + (2d) 2 



so that 



E(l) = -n 
V ; 2 




(5) 
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2d 2 
'b 



erfc 




Thus, in order to produce an approximate analytical ex- 
pression for the energy per electron E(v) in the LLL it 
is sufficient to derive values for the coefficients ctk in Eq. 
([3]). The effective thickness d* can then be calculated 
using Eq. ^. Of course, this approach cannot capture 
the small cusps in the energy associated with fractional 
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quantum Hall states, wh ich at finite temperature produce 
weak local maxima in rf** 6 * 7 * 30 l 

Our approach to estimating the coefficients at is as fol- 
lows. For a given value of d/ls and i/< 1 we calculate us- 
ing a numeric sum the Coulomb energy associated with a 
classical, triangular Wigner crystal of electrons interact- 
ing with the interaction law V(r). We also calculate the 
leading-order correction to this classical energy by using 
first-order perturbation theory for each electron in the 
slowly-varying Coulomb potential created by its neigh- 
bors. This potential is expanded to lowest order in the 
distance r from the potential energy minimum, and we 
use the LLL wavefunction ip(r) = exp[— r 2 /4Z^]/y / 27rIJ. 
The resulting estimate for E(y) is calculated for a par- 
ticular range of small filling factors, < v < v c , with 
v c -C 1, at which positional correlations are strong and 
the energy of the Wigner crystal state is a reasonable 
approximation for the total energy. Our choice of v c is 
discussed more fully below. Once E{y) is known in the 
range < v < v c , the function E(y) is then fit to the 
form of Eq. (|3} using a polynomial best fit with a^-o 
and ak>7 set to zero. The energy E(l) is given analyt- 
ically by Eq. mty. From the coefficients otk the inverse 
capacitance d*Jd is evaluated using Eq. The results 
of this procedure are shown as the black solid lines in 
Figs. [2} 

One reasonable estimate for the value of v c is v c = 0.2, 
since < v < 0.2 corresponds roughly to the Wigner 
crystal regime in the unscreened 2DEG 32 . While the 
screened 2DEG considered here may not remain in a truly 
solid phase up until v = 0.2, one can nonetheless expect 
that local positional correlations persist strongly, so that 
the WC state remains an accurate approximation to the 
total energy per electron. Fig. [2] uses v c — 0.2, but we 
find that our numerical results for d* are not noticeably 
different even if v c is made as small as 0.05. 



III. ANALYTICAL THEORY FOR TWO 
LIMITING CASES 

The general procedure described above allows one to 
numerically calculate d* within the LLL at arbitrary val- 
ues of v and d/ls- We now show that in the asymp- 
totic limits nd 2 ^ 1 and nd 2 <C 1, d*(v) can readily 
be described analytically. For the case nd 2 ^> 1, where 
the average distance between electrons is much smaller 
than d, electron-electron interactions are not significantly 
screened by the metal electrode. In this limit one can 
think that the system consists of a normal, unscreened 
2DEG with a uniform neutralizing background that is 
displaced by a distance d. By Eq. 0, the corresponding 
effective thickness can be written 



d*_ 
~d 



= 1 



rl 2 



2e 2 dd, 2 [VE ^ V)] 
when d/ls S> I and v, 1 



(7) 



where E-po is the FO expression^ for the energy per 
electron of the unscreened 2DEG with a coplanar neu- 
tralizing background. A form of Eq. Q was first 
proposed by Efros^ and used to explain the observed 
negative compressibility of 2DEGs in semiconductor 
heterostructureiP^. As we show below, however, Eq. ^ 
becomes incorrect if pushed to the limit where corrections 
to the geometric capacitance are large. 

We note that if one applies our method outlined above 
for calculating the coefficients ctk and the energy E(l) 
to the case of the unscreened 2DEG, one arrives at an 
expression for E(v) that is everywhere equal to E-po{v) 
to within 3.5%. Our semiclassical approach for estimat- 
ing E(u) at small v also reproduces the Hartree-Fock 
expressiorP^HUl U p ^ a term of order v 5 / 2 . Eq. ^ is plot- 
ted for a few different values of d/ls as the red dashed 
lines in Fig. [3j and indeed it matches our result in the 
appropriate limits d/ls 3> I and v, 1 — v ^> (ls/d) 2 . 

It should be noted that if Eq. (7| is naively pushed 
to the limit v <C (Is/d) 2 , it yields a negative capaci- 
tance. Such negative capacitance is forbidden by ther- 
modynamic stability arguments^, and this failure of the 
FO expression comes because at very small v the inter- 
action between electrons is screened by the metal gate^. 

Instead, in the limit nd 2 <C I the ground state can 
be described as a triangular Wigner crystal of electron- 
image charge dipoles, as shown schematically in Fig. [T] 
with interaction law V(r) ~ 2e 2 d 2 /er 3 . Treating these 
electrons classically and summing their interactions, as 
described above, gives an energy per electron E{v) ~ 
0.564(e 2 c? 2 /^ 3 3 )^ 3 / 2 . Adding the leading-order correction 
to this expression from first-order perturbation theory 
gives 



E(i 



fd 2 

l B 



(0.5641/ 3 / 2 + 0.429z/ >/z ) 



5/2n 



(8) 



when v <C min {l, /^/d 2 } 



This expression can be used to determine the coefficients 
a-k for 3 < fc < 7 by matching orders of v to Eq. p]), 
which gives: 



(9) 



v > l 2 B /d 2 , 



a 3 = 
a 4 = -E(l) 
a 5 = 0.564e 2 d 2 /2| 
a 6 = -2E(1) 
a 7 = lMe 2 d 2 /l 3 B 

when djls < 1 or f, 1 - 1/ < 1. 

Equations ([2| , (|3| , and (JsJ) can be combined to produce 
an analytical expression for d* . This is plotted as the 
blue dash-dotted lines in Fig. [3j (Due to its cumbersome 
algebraic form, we do not write this expression explicitly 
here.) For cases where d/ls < 1/2, d* is well described 
by Eqs. and ^ for the entire range < v < 1. 

For larger d/ls, this description is accurate only when v 
is very close to zero or unity. In the limit v — > 0, one gets 
d* jd ~ I.06(d//s)v / z/, which is equivalent to the classical 
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result derived in Ref. [T] Unlike in Ref. [TJ however, this 
result is correct even for 2DEGs with large Bohr radius, 
due to the confining effects of the transverse magnetic 
field. 
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FIG. 3. (Color online) The effective capacitor thickness, in 
units of d, as a function of filling factor in the LLL for various 
values of d/ls- The result of Eq. |7| is plotted as the red 
dashed lines and the blue dash-dotted lines correspond to the 
coefficients given by Eq. Q9B . The solid black lines are the 
result of a general numerical fitting procedure that determines 
the values of the coefficients at at each d/lB- 



Overall, our primary result shown in Figs. [2] and [3] is 
striking: it suggests that in the presence of a transverse 
magnetic field the capacitance can be arbitrarily large, 
limited only by temperature or the presence of a disorder 
potential. This is particularly surprising for graphene, 
where in the absence of magnetic field the capacitance 
becomes vanishingly gma lJ2Z128136| at ^ — > 0. This enor- 
mous enhancement of the capacitance in a magnetic field 
is a sister effect to Wigner crystallization in graphene, 
which is not possible in the absence of a magnetic field 
due to graphene's linear spectrum^, but which in the 
presence of a field becomes possible^ at small v. 



characteristic disorder impurity concentration Snd satis- 
fies Srid <C 1/d 2 . For the very thin devices examined in 
Ref. [231 where d ~ 5 nm, this implies Sn d «4x 10 12 

2 

cm . 

For devices with such small disorder, our predictions 
for d* should be correct up to some particular small value 
of v = is m , at which n becomes similar in magnitude to 
Srid, and d* rises sharply. This behavior is illustrated 
schematically in Fig. H (which resembles figure 2 of Ref. 
6) . A previous worH^I has shown that the minimum in 
the electronic compressibility for the unscreened 2DEG 
occurs atn« OASn^, which implies 



v m w 2.4(,5n d ); 



(10) 



Thus, in the window of filling factors u m < v < 1 — 
u m , the effective thickness is small and well-described by 
the theories presented here, while at smaller or larger v 
within the LLL d* is large. If the disorder is large enough 
(or the magnetic field is small enough) that (5rid)/B > 
3 x 10 10 cm _2 /T, this window closes completely, and the 
capacitance in the LLL is everywhere determined by the 
disorder. 
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FIG. 4. Schematic illustration of the effect of disorder on the 
measured inverse capacitance. The dashed line shows d* /d 
in the absence of disorder, and the solid line shows d* /d in 
the presence of a finite but small disorder potential, for which 
u m < 1. 



IV. DISORDER EFFECTS 



Of course, our analysis throughout this paper has ig- 
nored the effects of disorder, which truncate the diver- 
gence of the capacitance by destroying positional corre- 
lations. The presence of disorder modulates the density 
of the correlated electron liquid with some characteristic 
amplitude Srid- At small enough average electron den- 
sity that n < 5rid, pores open up in the 2DEG and elec- 
tric field lines starting at the metal gate leak through 
these pores. As a result, the electronic compressibility 
rapidly becomes small and d* grows sharply at small v 
and small 1 — v, as observed in Refs. [6| and [7] and illus- 
trated schematically in Fig. [3] Thus, observation of the 
large capacitance enhancement predicted here requires 
devices that are thin enough and clean enough that the 



V. CAPACITOR OF TWO 2D GASES 

So far we have focused on the case of a single 2DEG 
adjacent to a metal electrode. In this section we discuss 
capacitors made from two parallel 2DEGs, and show that 
similar capacitance enhancement effects can appear in a 
strong magnetic field. As a trivial example, one can no- 
tice that if one of the 2DEGs is sparse while the other 
has a much larger electron concentration, then this sec- 
ond 2DEG is equivalent to a metal electrode in terms of 
its electronic screening properties, and the theory out- 
lined above is applicable. On the other hand, one can 
consider a capacitor made from a sparse 2DEG and a 
parallel, 2D hole gas with equal hole concentration n (as 
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in Refs. 22-24). In such devices, for sufficiently small d 
(and in the LLL), electrons in one layer bind strongly to 
holes in the opposite layer, forming indirect excitons that 
interact as dipoleiP^. 

The electronic state of such a system was described by 
Yoshioka and MacDonalo® , who noticed that a sparse 
electron-hole bilayer system with total filling factor zero 
is equivalent to a bilayer electron system with total fill- 
ing factor one. That is, a system in which one layer con- 
tains electrons occupying a partially-filled LLL, v < 1, 
and the other layer contains an equal number of holes 
is equivalent to a system in which one layer has fill- 
ing factor v < 1 of electrons and the other layer has 
filling factor 1 — v < 1 of electrons. These authors 
showed that when the distance between the two layers 
is small enough that d/ls <C 1, the electrons in the 
two layers together effectively occupy the incompressible 
ground state wavefunctiorP^Vl/i i i, even though the con- 
stituent electrons are spread between two layers. This re- 
sult h as sinc e been confirmed by more detailed numerical 
studieJ^SBH. The corresponding pair distribution func- 
tion between electrons is therefore identical to that of 
the v = 1 incompressible state, Eq. GJ), so that the in- 
teraction part of the total energy per electron can be 
calculated as E- m t(v) = \n j d 2 rV(r)g(r), as in Eq. 
with V(r) = 2(e 2 /er - e 2 /eVr 2 + d?) a nd n = vl 2itij. 
Taking derivatives of this result for i<; int l 40 l 42 l 44 l 45 l gives, 
according to Eq. 



d* 
~d 



2~d 



1 — exp 



21% 



erfc 



lV2h 



In the limit d/ls <C 1, this expression becomes 



d* 

~d 



7T d 

81b' 



(11) 



(12) 



One can notice that this result for d*/d is similar in 
magnitude to the results presented in Sec. |III| for the mid- 
dle of the LLL; in both cases d*/d oc d/ls at small d/ls- 
In other words, capacitor devices comprised of two elec- 
tron or hole gases also exhibit capacitance that is greatly 
enhanced over the geometric value when d/ls *C 1 and 
the LLL is partially filled. As mentioned above, this 
result is unchanged if one considers electrons/ holes in 
graphene rather than in a conventional semiconductor 
quantum well where electrons have a quadratic disper- 



sion relation, due to the exact equivalence between the 
wavefunctions of these two systems in the LLL. 

Unlike the results shown in Fig. [2] Eq. ( 12 1 describes a 
constant, zn-independent capacitance in the LLL. On the 
other hand, when d/ls is of order unity, collective excita- 
tions about the ^1,1,1 state (electron correlations) lead to 
"softening" of the incompressible state and prod uces in- 
dependent corrections-^ to the result of Eq. (111. These 
corrections become larger with increasing d/ls, and at 
d/ls 3> 1 electr ons and holes in opposite layers decouple 
from each othei^^ 4 ^. A careful calculation of the ca- 
pacitance for double-2DEG capacitors at such d/ls > 1 
is left as the subject for a later work. 

Finally, we note that large capacitance enhancement 
may also be seen in situations where two parallel 
graphene layers with separation d are electrically con- 
nected (have the same electro-chemical potential) and 
are placed in a perpendicular displacement field D cre- 
ated by t wo dis tant metallic gates that are parallel to the 
graphene d 23 * 24 l In this case, the system of two layers is 
polarized along D. In a strong magnetic field this polar- 
ization takes the form of discrete, parallel, mutually re- 
pulsive dipoles, the charges of which are confined within 
disks of radius Ib within the graphene planes. Using our 
theory above one can find the 2D concentration of dipoles 
n as a function of the "applied voltage" V = \D\d, as well 
as the effective "capacitance per unit area" C — edn/dV, 
or, equivalently, the corresponding ratio d* j d. This ratio 
is given by exactly the same expression as for a real ca- 
pacitor made from two layers of graphene and charged by 
a real voltage (electro-chemical potential difference), as 
we described above. In the case of electrically conn ected 
graphene layers it is possible to deal with very smalP^] 
d and, therefore, potentially arrive at a giant capacitance 
enhancement. Substantial enhancement was already seen 
in Ref . E2 
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